How to measure the wavefunction absolute squared of a moving particle by using 

mirrors 
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We consider a slow particle with wavefunction ipt{x), moving freely in some direction. A mirror is 
briefly switched on around a time T and its position is scanned. It is shown that the measured reflec- 
tion probability then allows the determination of |V'T(a;)| 2 . Experimentally available atomic mirrors 
should make this method applicable to the center-of-mass wavefunction of atoms with velocities in 
the cm/s range. 



o 
o 

S3 



> : 

00 . 

ON ■ 

o ■ 
r» : 
o . 
\o ■ 
o ■ 

^ : 

9 L,: 
-i— > ■ 

G ■ 

s : 
cr 
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Measuring the absolute value of a wavefunction is part 
of the more general quest to reconstruct a quantum state 
from a series of repeated measurements on a system 

0, a 1 1 a a a a a m mmmmm- ^ ys - 

ically, the absolute square, |-0t(:r)| 2 , of the wavefunction 
of a particle gives the probability density for finding the 
particle at the position x, at time t. Hence, to determine 
it experimentally for an unknown wavefunction one just 
would have to perform a series of repeated position mea- 
surements. In general, the difficulty with this procedure 
is the required precision since it must be much better 
than the spatial extent of the wavefunction. An addi- 
tional problem occurs for slow laser cooled atoms since 
their wavefunction can be influenced by a measurement. 
Here we are mainly interested in the center-of-mass (CM) 
motion of such slow laser-cooled atoms with wavefunction 
widths in the region of a few micrometers and velocities 
in the cm/s range. Such states can be obtained by letting 
atoms escape from a trap (16|. 

An elegant position measurement technique [13] , which 
would have the necessary resolution, places the atom in 
a strong gradient of a magnetic or light field. In this 
field the internal energy levels vary rapidly with posi- 
tion. By sending a light beam at one of the atomic reso- 
nance frequencies one then detects the atoms which are 
located where the resonance occurs. In this way, sub- 
wavelength resolution can be achieved. However, when 
the atom moves through the field its CM wavefunction 
will, except under special circumstances, be affected by 
the interaction and will not remain the original free CM 
wavefunction. Hence one might not obtain ^(a;)! 2 but 
some other distribution. 

In this paper we propose an alternative method for de- 
termining |-0 t (a;)| 2 for very slow free particles, in partic- 
ular for laser cooled atoms. This method uses reflection 
by a mirror which is switched on for a brief time period r 
and whose position is scanned. The reflection probabil- 
ity which can be determined sufficiently far away from 
the mirror by fluorescence techniques, will be shown to 
allow the determination of | -t/>* (as ) | 2 . 

A suitable mirror is an evanescent wave (EW) mirror 
[l9L I20I l2~fl ] , depicted in fig. Q] A laser beam enters 
a prism and is totally reflected at the surface, creating 
an evanescent wave outside the prism. Due to the in- 
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FIG. 1: Schematic scetch of an evanescent wave mirror. The 
atoms are reflected by the potential created by the evanescent 
wave which occurs on the surface of the prism. 



tensity gradient an atom positioned there experiences a 
steep exponential potential which is repulsive if the laser 
frequency is above the atomic resonance frequency (blue 
detuning). For a strong laser, a slow atom will be re- 
flected without reaching the surface. Such a mirror is 
easily switched on and off by an acousto-optic modulator. 
Switching times of .5 (is or less can be implemented. Re- 
flection takes place within a region in the order of 0.1/im. 
For our purposes the experimental setup could be similar 
to the temporal two-slit experiment of Ref . [l6| . 

To simplify the analysis we replace the exponential po- 
tential of the mirror by a high square potential barrier. 
Letting its height go to infinity and its thickness go to 
zero the mirror is modeled by a reflecting plane. Possi- 
ble adsorption on the prism of the mirror is disregarded. 
This modeling simplifies the mathematics but keeps the 
essential physics. It will be shown that for small r the 
reflection probability scales as y/r, multiplied by a factor 
depending on the absolute square of the wavefunction at 
the position of the mirror. It is also shown that for larger 
r there is a region of linear dependence. 

First we consider the one-dimensional case, with a free 
wavefunction coming from the left and traveling to the 
right (see fig.H]), which is denoted by 'tpt(x) if there is no 
mirror present. When the mirror is active from T~t/2 to 
T+t/2 the wavefunction evolves freely until time T— r/2, 
experiences the mirror from T — r/2 to T + t/2 and then 
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FIG. 2: Wave packet ipt(x) coming in from the left (solid line). 
Mirror (dashed) at xm switched on at T — r/2 and switched 
off at T + r/2. Partially reflected wave (dotted). Measured 
reflection probability allows the determination of |V>t(:em)| 2 . 



again evolves freely. The mirror action causes partial 
reflection, and the reflected part can be calculated either 
by an eigenfunction expansion, or with Green's functions 
as in the following. 

The mirror is first positioned at xm = 0. The free 
time development operator (7 F (i,0) is given in terms of 
the free Green's function G F (x) by 

(x\U F (t,0)\x>) 

= Gf(x - x') = y/^j2^me- mix - x ' )2/2iht . (1) 

During the time the mirror is active the time development 
operator U^(t, 0) on the right (+) and left (— ) half axis, 
respectively, is given in terms of the Green's function 
G^Oby 

{xlT^faOW) =Gf(x.x') 

= 6±(x)(Gf(x - x') - G F (x + x'))9 ± (x') 

where 0±(x) = 9(±x) and 9{x) is the Heaviside function. 
Let (f)f denote the wavefunction on the right and left 
half axis, respectively, when the mirror is active. With 
"0t-t/2 the (free) wavefunction on the whole axis at time 
T - r/2, one has 4^ +t/2 = U ± (t, 0)(9±ip T _ r/2 ), i.e. 

= 0±(x){G F T * (0±^_ r/2 ))(x) 

-9 ± (x)(G F T *Tc(6 ± ij T _ r/2 ))(x), (2) 

where * denotes convolution and 7r reflection: x — > —x. 
The complete wavefunction at the switch-off time is 
4>t+t/2( x ) — ( t ) T+T/2^ x ) + < ^t+t/2( 2 ')' wmcn then evolves 
freely at later times. The reflected part is given by the 
Fourier transform 4>T+T/2{k) for fc < and the reflec- 
tion probability, denoted by iV re fl(r, T; Xm), is given by 
its norm squared, i.e. for Xm = 



iV refl (T,T;0) 



dfc|0 T+r/2 (fc)| : 



(3) 



Fourier transforming eq. ^ one obtains after some cal- 
culation 

4>T+r,2{k) = ~(^T-r/ 2 (*0 " £r_r/ 2 (-fc)) ^ 
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where & denotes Cauchy's principal value and 
a = tvr /2m . 



(4) 



(5) 



For t — ► the right hand side of eq. ((4| has to converge 
to ipT(k) which shows that there must be compensations 
between the first and second term. As a consequence, the 
behavior for small r (i.e. small a) is not easily extracted 
from this equation. We therefore use the skew-reciprocity 
of the Hilbert transform L 22j to write 



l(jf T - T /2(k) - ^T-r/2(~k))e- iak2 



i>T-r/2(k)e 
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dfc 



k' - k k' + k 

„ ipT-r/ 2(k") iak „2 
k" 



(6) 



and then use the anti-symmetry in k' of the square 
bracket in eq. (j4]) to obtain, with the substitution k = 

fc'/VS and $r_ T/2 (fc) = e iQfe2 /V T (fc), 

0T + ./ 2 (fc) = Mk)e~ iak2/2 



+ 
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dfc 
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dfc" 2k- 



ak" 2 



i iak " 2 / 2 ib T {k"). (7) 



The last, now non-singular, integrand can be expanded 
in \fu, i.e. y/ak <C 1 for typical fc values of ?/>, which in 
lowest order yields for eq. (J7J) 



T+T/2 (fc) = Mk)e~ iak2/2 

In r<x> 

+ V^-Vt(0)^ / dn 



1 



- 1 



k — \fak 



(8) 



where / &k' iprik") — V2ttiPt(0) has been used. For 
fc < one has tpT(k) = since the free particle is com- 
ing in from the left. The reflected wavefunction has a 
yfa. in front and from this one might, erroneously, expect 
that the reflection probability goes as a, i.e. as r. How- 
ever, calculating the reflection probability with Eqs. ([3]) 
and ([5]). while using the symmetry of the appearing in- 
tegrands and employing footnote [22j], the substitution 
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k — > k/y/a gives an additional X/y/ot and one obtains, 
with the mirror located at im = x, 



(9) 



N Ie& (T, T; x M = x) = 2\ \ifa(x)\ 



in lowest order in y/r, with the mirror active from 
T - r/2 to T + r/2. Since / dx |?/>t| 2 = 1, eq. © im- 
plies 

\Mx)\ 2 = N icR (t,T;x m = x)/ { dx' N IeR (T,T;x'). 



(10) 

Going to higher orders in the expansion of the last inte- 
gral in eq. ([7]) one obtains after some calculation, again 
using the symmetry of the appearing integrands and em- 
ploying footnote |22j |. 



N ie s_(T,T;xM = x) =2 

tlT 



— Re (^r(x)^(a?)) + 0(7 
Dm 



(11) 



Thus eq. ((9]) is expected to be a good approximation if r 
satisfies r -C 6m|-0T(xM)|/^|'0T'(xM)|. This shows that if 
the curvature of the wavefunction is large one needs short 
switch-on durations of the mirror for eq. ([9]) to be a good 
approximation. Equation (|10[) , on the other hand, turns 
out to be valid for a much wider range of mirror pulses 
r. 

It will now be shown that there is a linear regime for 
larger r and that eq. (TT0|) also holds in this regime. By 
means of Eqs. ([3]) and J7]) one can write AT re fl in the form 



7V rcfl (T,T;0)= I dk 1 dk 2 Mh)I(ki,k 2 )Mk2) (12) 
where a lengthy calculation gives for the integral kernel 

I(k u k 2 ) 



x < e 



{' 



2tt 

-io(fe^-fc?)/2 



kl) 



[fcierf(-s/iafci) + k 2 ert(y/— iak 2 )] 



- e iQ ( fc 2-fc?)/2 [ fc2 e rf(Vkwfc 2 ) + fci exf(\/=iaA;i)] }. (13) 

Expanding this for small a, i.e. -y/afc <C 1 for typical k 
values of ip, one recovers Eqs. © and fTT]). On the other 
hand, when -^afc 3> 1 one can use the asymptotic expan- 
sion of the error function, erfz = 1 + u^ 1 / 2 e~ z2 {— X/z + 
l/2z 3 + (D(z~ 5 )}. With this one obtains in lowest order 

t,i , a /, , \ sina(fc| — fc?)/2 

/(fc 1 ,fc 2 ) = -(fc 1+ fc a ) a{ ^_ k 2y2 ■ (14) 

The contributions from the next order cancel and z~ 3 
yields an a -3 / 2 term. Provided afcoA/c is small, with hk 
the mean momentum and HA^ the momentum width, the 
last fraction in eq. (|14|) can be approximated by 1 and 




FIG. 3: Reflection probability vs. pulse duration r for a slow 
(vq = l.lcm/s) initially minimal Gaussian wave packet for 
Cs atoms prepared at time t = with mean position xo = 
— 0.66 /im, initial width A x = 0.1 /im. Mirror switched on 
at T = 60 /is. Solid line: numerical result. Dashed line: 
lowest order (eq. ©). Dashed-dotted line: up to second order. 
Dotted line: up to third order. 

FIG. 4: Wavefunction and T as in fig. EJ Solid line: \f T (x)\ 2 . 
Dashed-dotted line: iV rc fl (r, T; Xyi = x) /2y/hr firm for r = 
1 /is. Dotted: the same for r = 5 /is. Fat dots: r.h.s of 
eq. (|10[) for r = 1 /is. Dashed: r.h.s of eq. (|10p for r = 5 /is. 



k\ and &2 by fco- Then, in this order, 7V rc fl(r, T; xm = 
i) = 2|^ T (a;)| 2 fco « = |V't(^)| 2 wo''" j which yields eq. (fTU|> . 
Under the same assumptions the next nonvanishing cor- 
rection is small. Hence in this regime the reflection prob- 
ability is linear in r and eq. (|10[) is a good approximation 
[23l | . For the wavefunction considered in fig. [5] one has 
y/ak = 25 3> 1 and afco^fc < -8 which is sufficient in this 
case. 

In fig. [3] a numerically calculated reflection probability 
-/Vrefl(T, T; xm = 0) as a function of r is compared with 
various analytic approximations. For the chosen example 
of a Gaussian wavefunction of a Cs atom with velocity 
vq = 1.1 cm/s it is seen that eq. (|9|) agrees well with 
the exact result for r < 1.5 /is. For higher velocities 
shorter pulse durations r are needed in eq. ^ . If these 
are too short for available mirrors one may either fit the 
higher order result in eq. (jlip to experimentally obtained 
reflection probabilities for different pulse durations r or, 
better, use eq. (|10p in the regime linear in r. 

In fig. 2] we consider the same wavefunction as in fig. [3] 
The solid line is \iPt{x)\ 2 ■ This is compared with the nu- 
merically obtained reflection probability, iV re fi(r, T; im = 
x), for different mirror positions x, divided by 2^/ fvr j nm 
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FIG. 5: Slow (vo = 25 cm/s) initially minimal Gaussian wave 
packet for Cs atoms prepared at time t = with mean position 
xq = —12 mm, initial width A x = 1.6 fim. Mirror switched 
on at T = 50 ms. Plotted: |i/>t(:e)| 2 (solid line) and r.h.s of 
eq. (fTTI)) for r = 10 (is (dots) and r = 30 (is (dashed line). 

in view of eq. © . For r = 1 (is (dashed-dotted line) the 
agreement is good, but not quite so good for r = 5 (is 
(dotted line), as expected from fig. [3l Also in fig. [4] is 
plotted the right hand side of eq. (flU)) for t = 1 (is and 
r = 5 (is. For r = 1 (is (fat dots) it is practically in- 
distinguishable from 1^(^)1 2 and now also for r = 5 (is 
(dashed line) the agreement is excellent. The slight shift 
is due to the next order term in eq. (Til"]) . 

A test of eq. (1101) in the linear regime is displayed in 
fig. [5] for a slow (vq = 25 cm/s) wavefunction for which 
the lowest order expressions in Eqs. (J9j> and ([8} are valid 
only for unrealistically short pulse durations of less than 
0.01 (is. However, even for pulse durations of up to r = 
30 (is the right hand side of eq. (|10p agrees extremely well 
with \ipx(x)\ 2 . This is remarkable since during 30 (is the 
wavefunction moves by 7.5 /im, a third of its width at 
time T, and the agreement is attributed to interference. 

In the three-dimensional case the mirror is first as- 
sumed to be positioned at x m and parallel to the y — z 
plane. The wave packet is assumed to be incident from 



the negative x direction, i.e. ipt(k) is nonzero for k x > 
only. The analysis is analogous to the one-dimensional 
case and one obtains 

/ dy J dz \i/j T (x M ,y, z)\ 2 

= N le a(T,T;x M )/ J dx' N TeS _(T, T; x) (15) 

both in the square-root and linear regime. By this ex- 
perimental setup one therefore obtains the position space 
density for the wavefunction integrated over the mirror 
plane |24l |. By rotating the mirror by different angles 
around the point xm one can also obtain the integrated 
atomic density in these rotated planes. In addition one 
may also vary the point xm- This is similar to the situ- 
ation encountered in computer tomography, and in this 
way one can recover IV'tC^c) | 2 from the experimental data 
by standard techniques. The presence of a perpendicular 
gravitational field does not change the conclusions pro- 
vided the pulse duration is sufficiently short. Indeed, un- 
til the mirror is switched on the free evolution is modified 
by the gravitational potential, but during the short pulse 
time its effect is negligible and the reflection probability 
is similar as before. 

In summary, we have proposed to measure the absolute 
square of the wavefunction of a slowly moving particle by 
using reflection from a pulsed mirror. It has been shown 
that for short pulse duration r the reflection probabil- 
ity scales as \fr and for somewhat longer duration as 
r. By measuring the reflection probability and scanning 
the mirror position one can determine |-i/>i (cc) | 2 . With 
presently available atomic mirrors it should be possible 
to apply this method to the center-of-mass wavefunction 
of atoms with velocities in the cm/s range. This tech- 
nique provides a sub-wavelength resolution and has the 
advantage of avoiding, at least in principle, the problem 
of a possible disturbance of the wavefunction during the 
measurement. 
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